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Bipartite entanglement in the ground state of a chain of N quantum spins can be quantified either 
by computing pairwise concurrence or by dividing the chain into two complementary subsystems. 
In the latter case the smaller subsystem is usually a single spin or a block of adjacent spins and the 
entanglement differentiates between critical and non-critical regimes. Here we extend this approach 
by considering a more general setting: our smaller subsystem Sa consists of a comb of L spins, 
spaced p sites apart. Our results are thus not restricted to a simple 'area law', but contain non-local 
information, parameterized by the spacing p. For the XX model we calculate the von-Neumann 
entropy analytically when N —* oo and investigate its dependence on L and p. We find that an 
external magnetic field induces an unexpected length scale for entanglement in this case. 



PACS numbers: 03.65.Ud, 03.67.-a, 73.43.Nq, 75.10.Pq 



INTRODUCTION 

Quantum phase transitions at zero temperature cor- 
respond to a fundamental restructuring of a system's 
ground state. In quantum spin chains these transitions 
occur as an external parameter (e.g. a magnetic field) 
is varied [l|, and are manifested as a marked change in 
the decay of quantum correlations: algebraic at the crit- 
ical point and exponential away from it. The amount 
of entanglement present in the ground state is expected 
to depend significantly on whether the system is critical 
or not, since at a critical point all the constituent parts 
of the system must be non-locally correlated and thus 
entangled 11111111110. 

Unfortunately, it is not yet clear how to measure entan- 
glement in general l8|. At present, we only understand 
fully how to quantify bipartite entanglement |9|. It is 
natural, therefore, to try to extract as much information 
as possible about entanglement in this context. Ways of 
doing this have recently been the focus of considerable 
attention. One possibility is the following. From a chain 
of N spin-i particles select two, compute their reduced 
density matrix, and then obtain the associated concur- 
rence 10]. This is a function of the separation between 
the selected spins. Despite the fact that the concurrence 
is only a short-range measure (it vanishes if the spins are 
farther apart than next-nearest neighbors), this approach 
has been applied to detect phase transitions in a variety 
of situations 0j S ■ A second possibility is to measure 
the entanglement between a single spin and the rest of 
the chain. This has also been related to the presence of 
a critical point J3]. 

The problem with these methods, which involve a small 
number of spins, is that they do not take into account the 
fact that entanglement is shared between many parties, 
i.e. they provide little information regarding the non-local 
nature of entanglement. This deficiency is shared by an- 
other much-studied bipartite division of the spin chain, 
namely that between a block of L adjacent spins and the 
remaining iV — L [J,!^. In the limit N ^ oo, the entan- 



glement entropy has in this case been computed analyt- 
ically using the theory of Toeplitz determinants (for the 
XX model) |llLll2| . conformal field theory [l3j, and from 
averages over ensembles of random matrices l],4|. For 
one-dimensional chains it has been shown that as L — > oo 
the entropy tends to a constant value away from critical 
points, and that it diverges like InL at phase transitions. 
For critical d-dimensional spin-lattices this 'block' entan- 
glement has been pro ven to grow like L'^~^\nL under 
certain conditions J15l|. while for a gapped system one 
expects an area scaling law due to the finite correlation 
length jig . Such a direct relation between entanglement 
and area is known to hold in harmonic lattices [13 • All 
these results indicate that, at least for large blocks, the 
entanglement comes mostly from the boundary. 

Our purpose here is to introduce a new geometry in 
which to study bipartite entanglement in quantum spin 
chains. We divide the chain into two subsystems, Sa 
and Sb, as follows. Sa consists of L equally spaced spins, 
such that the spacing between the spins in this subsystem 
corresponds to p sites on the chain. Sb then contains the 
remaining N — L spins. (Obviously this only makes sense 
if A^ > {L — l)p.) Sa can be visualized as a comb with 
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FIG. 1: The comb division illustrated for three different values 
of the spacing p. In all cases one subsystem {Sa, denoted by 
black circles) contains L spins and the other {Sb, denoted by 
empty circles) contains the rest of the chain. We have written 
explicitly the label of the last of the L spins (the label of the 
first spin is 0). The first case, p = 1, corresponds to the well 
known 'block' division. 



L teeth. This geometry enables us to study non-local 
entanglement effects by varying the spacing p. 

Three possible divisions are illustrated in Fig.l. For 
p — 1 we recover the simple 'block' arrangement, which, 
as was noted above, has already been the subject of ex- 
tensive investigation. In this case the subsystems Sa and 
Sb are only 'in contact' near their common border, and 
as a result, for large values of L, the dependence of the 
entanglement on L is minimal. For p > 1 the entangle- 
ment between Sa and Sb is shared between many dif- 
ferent sites and it then grows linearly with L, to leading 
order as L ^ cxd. The logL term that dominates when 
p = 1 appears as a secondary correction when p > 1. As 
an example, we investigate how the comb entanglement 
in the XX model depends on the spacing p when L is 
fixed and show that this reveals the emergence of a new 
length-scale determined by the external magnetic field. 



BLOCK ENTROPY 

The XY spin chain in an external uniform magnetic 
field h has as its Hamiltonian 
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where {aj ,(j^ ,(7^) are the usual Pauli matrices. When 
7 = this system is called the XX model and when 
7 = 1 it is the Ising model. It is an integrable model [l3| 
which displays both critical and non-critical regimes. In 
the 7 — ft. phase diagram the lines ft, = ±1 and the seg- 
ment 7 = 0, |ft| < 1 are critical, while all other regions 
are non-critical. For simplicity, we will restrict our anal- 
ysis mainly to the XX model and will come back to the 
general case towards the end. 

Since the XX ground state is just a (non-entangled) 
ferromagnet for |ft| > 1, we will only consider its critical 
regime. A more convenient parameter in that case is the 
angle k defined by 



ft = cosfc, k e [0, tt]. 



(2) 



calculation of the expectation value of a product of any 
number of JW operators. 

The matrix Cn factorizes into a direct product, 



Cn = T[g] 
where T[g] is the matrix 

iT[g])jk = 9j-k, gn = 



1 
-1 



(4) 



27r 



2tt 



g{e)e-"^'de. (5) 



The function g{0) is called the symbol of the Toeplitz 
matrix T[g\. For the XX chain it is given by 
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1 li -k<e <k, 

— 1 otherwise. 



(6) 



Following the calculation presented in [ll| , the entropy 
of subsystem Sa is obtained as a contour integral in the 
complex plane: 



E = lim lim 

e^0+ 5^0+ ZTTl 



c{e,5) 



..dhlDAW 

e(l + e, A) — dX, (7) 



where Da{X) — dct{XI —TA[g]) involves the matrix rA[g], 
which is obtained from the original matrix T[g] by remov- 
ing the rows and columns that correspond to sites in Sb , 
and 

x + y fx + y\ x-y (x~y 

e{x,y) = -^log2 (^^j - ^-lo& [-^ 

(8) 
The contour of integration c(e, 6) approaches the inter- 
val [—1,1] as e and S tend to zero without enclosing the 
branch points of e(l -I- e, A). 

When Sa corresponds to L consecutive spins, i.e. in the 
block case, Ta is simply a block inside T and thus is also 
a Toeplitz matrix. This allowed Jin and Korepin |ll| 
to obtain the corresponding entropy by using a proved 
instance of the Fisher-Hartwig conjecture relating to the 
asymptotics of Toeplitz determinants |l2| . This states 
that in the limit of large blocks we have 



InD^(A) =coL + (3^liiL + 0{l), 



(9) 



We denote by \'^) the ground state of this system and 
introduce the Jordan- Wigner ( JW) transform at each site 
of the lattice. 
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Given any product of an odd number of JW operators, 
we can see from the symmetry of the Hamiltonian that its 
expectation value with respect to JVP) vanishes. Wick's 
theorem and the relation {^\mjmk\'i') = Sjk + i{CN)jk, 
where Cjv is called the correlation matrix, allow for the 
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(10) 



The coefficient f3{X) may be calculated by writing the 
symbol as 

A - g{9) = HX)t[3{9 - k)t-f3{9 + fc), (11) 

with t0{9) = e-^P'^"'-"^ and 



<^(A) = (A+1) 



A + 1 
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(12) 



and hence it is given by 



Its values are 



«^)-^'"(^ 



(13) 



As a consequence of the Fisher-Hartwig conjecture the 
entanglement as i — > cxd is 



E{L) =£iL + £2\nL + 0{l) 



(14) 



Because of the simplicity of the function g{6) all the rel- 
evant quantities can be evaluated. In particular, we have 



co(A) = -[/cln(A-l) + (7r-A:)ln(A + l)], (15) 

and by substituting this into Q we see that the leading 
order term £i actually vanishes because 



e(l,l) = e(l,-l) = 0. 



(16) 



Therefore the dependence upon L is only logarithmic, 
with a prefactor £2 which is obtained from 0. It turns 
out that when p > \ the leading order contribution does 
not vanish. The next section is devoted to its computa- 
tion. 



COMB ENTROPY 

The choice we propose for the subsystem Sa also leads 
to a Toeplitz structure for the matrix Ta- We single out 
those spins whose label is a multiple of an integer p, and 
thus we have 

{TA),k = 9pj-pk = ^ fj g{e)e-^^^^-''^'>de. (17) 

This is not yet in the Toeplitz form. We must first find 
a function gp{9) such that 

g(a)e-'P""da = / gp(a)e-""da, (18) 

and this will be the symbol of Ta- Multiplying H18|l by 
gjnS .^yj^]-^ < < 27r, summing over n and using the 
Poisson summation formula we arrive at 
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(19) 



so the value of gp at the point 9 is obtained as the average 
value of g over the p vertices of a regular polygon. 

This average is not hard to calculate. It is easy to see 
that for each k the function gp{9) is piecewise constant 
and even, with jumps at the critical points ±0* given by 



[0,7r) 9 6** = min{a,27r-a}, a =pk mod 2tt. (20) 
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(21) 



where the brackets |-] denote the integer part and 

s = signja — tt}. (22) 

The entanglement will now depend on the spacing: 

E{L;p)^£i{p)L + £2{p)liiL + 0il) (23) 

as L ^ 00. To calculate it to leading order we only need 
the integral 
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(24) 



and this leads to 

£,ip) = TT-^ [re(l, 5p(0)) + (vr - r )e(l, gpiir))] , (25) 

which is our main result. We calculate the logarithmic 
correction £2{p) in the next section. Note that the bound 
< fi(j') < 1 is explicitly respected. 

If we restrict our subsystem to be a single spin, i.e. L = 
1, it is easy to calculate the entanglement because the cor- 
relation matrix has only one element, {2tt)~^ L^ g{9)d9. 
We get simply 
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(26) 



On the other hand, from equation H21(l we obtain that as 

p — > 00 



and 
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(27a) 



(27b) 



where both qi{p) and (72 (p) vanish like p~^ . If now we 
insert (|27|l into the general expression H25|l , then we find 
that the total entanglement in the limit of large spacing 
converges, as expected, to a combination of single-spin 
contributions. 



fi(p) 



P 



(28) 



where a is a constant. It is interesting to note that the 
rate of convergence is rather slow, indicating a long-range 
dependence of the entanglement on the spacing. This is 
in contrast with the behavior of simpler quantities like 
the pairwise concurrence, for example, which vanishes if 
the spins are more than two sites apart. 
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FIG. 2: (color online) E\ as a function of k for various values of 
the spacing p. The curve is always piecewise linear and weakly 
converges to E\ . Notice that for k — 7r/2 the entanglement is 
maximal whenever p is even. 
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FIG. 3: (color online) £i as a function of the spacing p for 
various values of the magnetic field h = cos k. For k = tt/^ 
there exists a typical length scale for the spacing: the entan- 
glement is maximal for p — n£, and its value at these points 
is equal to the limit Ei. 



We next consider the case when p 
(72(0) and (?2(7i') are given by 



52(0)-^(.g(0) + 5W) 







2. The values 



(29a) 



and 



ff2W = -{g{n/2)+g{3n/2))=g{7r/2). 



(29b) 



The latter is either 1 or —1, and thus makes no contribu- 
tion to entanglement. The critical angle is just 6* — 2k 
and hence we have fi(2) = 2Lk/'K. The linear depen- 
dence on k can be seen in Fig. 2, where we plot the en- 
tanglement for various values of the spacing p. In the ab- 
sence of any external magnetic field, i.e. for k = 7r/2, the 
entanglement attains its maximum possible value when- 
ever p is even. For larger values oip the function is always 
piecewise linear, eventually converging to Ei. 

It is important to observe that the expression (|25|l for 
the entanglement is continuous when we consider p as 
a real number, despite the discontinuities that appear 
in 1)21(1 . The function 5p(0) is discontinuous whenever 
pk — 2mr, but at those points 9* vanishes and hence 
£i{p) remains unaffected. On the other hand, £i(j>) is 
also oblivious to the jumps in the function gpij^), which 
occur at pk = (2n + l)7r (due to the variable s), be- 
cause then we have 6* = tt. Its derivative, on the other 
hand, is discontinuous: at the special points pk = rnr the 
entanglement has a local maximum with a cusp form. 
Remarkably, for a fixed k its values at these maxima are 
all the same (i.e. they do not depend on n) and are equal 
to the large-p limiting value. 



In Fig. 3 we plot £i{p) as a function of the spacing for 
different values of the magnetic field h = cosfc. Since, 
of course, only integer values of p may be realized in the 
actual chain, we take k = ir/i, where £ is an integer. We 
see that this leads to the appearance of an unexpected 
length scale for the entanglement: it attains its maximal 
value whenever the spacing is a multiple of £. The exis- 
tence of such a length scale appears to be a fundamental 
property of quantum spin chains in magnetic fields. 

Recently the formalism of Toeplitz determinants has 
been used to compute the 'block' entanglement of the 
more general XY model J19j, with finite anisotropy pa- 
rameter 7^0. The main difference with respect to the 
case 7 = is the form of the function g{9), which is no 
longer piecewise constant and becomes complex. Calcu- 
lating the average (|19|l then becomes much less simple, 
but it is still possible to employ the present approach to 
obtain the entanglement for arbitrary values of the spac- 
ing p. A special case for which explicit calculations are 
possible is the Ising model without magnetic field, ob- 
tained from QJ by setting 7=1 and h = 0. In this 
case we have g{9) — e'^ and thus gp{9) — 0, leading to 
the result that EisingiL;p) — L for any value of p. This 
reflects the fact that the zero temperature ground state 
of this model is maximally entangled [3| . 



LOGARITHMIC CORRECTION 

In order to obtain the logarithmic correction for the 
XX model we need to decompose the symbol (fTI^ in the 
form 



£j^{mr/k) ^Ei. 



(30) 



A - 5j>W = H^)U9 - 9*)t-pi9 + r), (31) 
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FIG. 4: Left: £2 as a function of k for two values of the 
spacing. Right: £2 as a function of p for two values of the 
magnetic field. 



where 



'/'(A) = (A - gp(7r)) 



A-ffp(0)y 



(32) 



is independent of 6 and the discontinuities are accounted 
for by 



The function /3(A) is given by 



/3(A) = -iln^^ 



2Tii VA-gp(O)/' 



(33) 



(34) 



and the coefhcient in the logarithmic correction to the 
entanglement is 



S,{p)^-<fe{l,X)(3{X)^dX. 
in / dA 



(35) 



Since for p > 1 both |(?p(0)| and |gp(7r)| are smaller 
than unity, we can choose as our contour of integration 
the unit circle. Using power series expansions we arrive 
at 
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(37) 



Since both gp{0) and gpiir) are piecewise constant as 
functions of k, the same is true for £2- In the left panel 
of Fig. 4 we see that the number of jumps in this function 
grows as p increases. On the other hand, for a given value 
of the magnetic field £2 decays rapidly as p increases, 
and has discontinuities at p = ni when k — n/i (see the 
right panel of Fig. 4). Notice that the vanishing of the 
logarithmic term in the entanglement is consistent with 
the fact that E -^ LEi as p — > cx) (cf. IJ^SI)- 



CONCLUSIONS 

In summary, we have extended the bipartite approach 
to entanglement in one-dimensional critical spin chains 
by introducing a new partition of the chain, the comb 
partition, which allows us to go beyond the simple 'block' 
picture and investigate non-local correlations analyti- 
cally. The organizing subsystem consists of L spins sep- 
arated by p sites, and we have found that as p ^ 00 its 
entanglement with the rest of the chain reduces to the 
sum of the individual contributions of its elements, al- 
though with a slow convergence rate that indicates the 
existence of long-range correlations. We have also found 
that the presence of a magnetic field induces a typical 
length scale for entanglement. It would be interesting 
the see if this length scale is present in other statistical 
properties of critical spin chains. Our results regarding a 
generalized version of the Emptiness Formation Probabil- 
ity, which has recently been c omp uted for the XY model 
using Toeplitz determinants 2(J|, will appear elsewhere 
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